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Partition Function Zeros in Non-compact QED 

A. Ali Khan and I. Barbour 
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The lowest zeros of the lattice partition function for non-compact QED are found in the complex fermion mass 
plane on 6*, 8* and 10* lattices at intermediate values of the coupling. The scaling of the low lying zeros with 
lattice size is analysed. 
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1. Introduction 

We present preliminary results for the be- 
haviour of the lowest zeros in the complex fermion 
mass plane of the lattice partition function of 4- 
dimensional non-compact QED with four flavours 
of staggered fermions. Results on G'*, S'* and 10'* 
lattices are presented for values of the inverse cou- 
pling /3 in a region where there is a possible tri- 
critical point in the (zero temperature) infinite 
volume limit. 

The work presented here is the first part of an 
ongoing study of the scaling behaviour of the ze- 
ros as a function of increasing lattice volume. 

2. Method 

An ensemble of configurations is generated us- 
ing a Hybrid Monte Carlo algorithm, choosing 
the molecular dynamics parameters to keep the 
length of each molecular dynamics trajectory ~ 
0.8 and the acceptance rate between 0.65 and 0.8. 

Measurements are made on configurations sep- 
arated by approximately 3 or 4 units of molecular 
dynamics time. 

The lattice partition function is defined as 
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where mo is the update mass, M is the fermion 
matrix at zero fermion mass, m, and 5*^ is the 
standard non-compact gauge field action. The 
eigenvalues of M are found using the Lanczos 
algorithm [1] to tridiagonalise M. The tridiag- 
onal form is then diagonalized using Sturm se- 
quences. Since iM is hermitian, this method gives 
the eigenvalues without reorthogonalization and 



its consequent storage demands. Accurate deter- 
mination of all the eigenvalues enables accurate 
determination of the coefficients of the character- 
istic polynomial (each weighted with respect to 
the determinant at the update mass) [2]. 

These coefficients, averaged over the ensemble, 
give the polynomial expansion in for the par- 
tition function 
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on an lattice. 

We choose the exponential parametrization of 
the coefficients because they vary within a large 
range of values (on a G'* lattice between about 
and e^^° , on a lO'* lattice between about and 
g4400-j ^j^j gj.Q^ linearly with the volume. 

The zeros of this polynomial are the Lee- Yang 
zeros in the complex mass plane. 

Since the eigenvalues of M appear as complex 
conjugate imaginary pairs, Z is a polynomial in 
and its coefficients are necessarily positive. 
The fact that the coefficients vary grossly in mag- 
nitude gives rise to technical difficulties: 

a) in their determination 

b) in finding the roots of the polynomial. 
These difficulties have been overcome and the 
methods adopted will be described in a subse- 
quent publication. 

3. Results 

For non-compact QED with staggered fermions 
on a G'* lattice for 0.16 < (3 < 0.22, the zeros near 
the physical mass region (m > and real) are 
imaginary. This suggests that any phase tran- 
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sition will only occur at zero fermion mass, as 
expected. Note that this signal differs from that 
in compact QED [3], where the lowest zeros at a 
coupling of /3 = 0.885 are found to be complex 
with non-zero real part. 

Errors on the zeros are not shown since they 
are difficult to estimate in a rigorous manner. The 
statistical errors can be estimated via a bootstrap 
technique. This gives an error of about n% on the 
n'th zero (with n=l corresponding to the lowest 
zero). A related systematic error will also arise 
from limited ensemble sampling. We estimate 
this error to be comparable with the statistical 
error. 

It is also relevant to note that the a„ are ap- 
proximating to a continuous function of n. The 
equation 
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(k integer), becomes, if the a„ become continuous 
in the infinite volume limit, for V large 
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and its roots are the same as those of Z . That 
is, using equation 3 with k = 1,2,3,..., we find 
the same roots near the physical region. This 
signal for the continuity of a„ grows stronger with 
increasing lattice volume. 

Figure 1 shows the lowest zeros plotted in se- 
quence. It is easily seen that in the infinite vol- 
ume limit and in the limit m ^ 0, the condensate 
is given by p(0), the density of zeros at the origin. 
For large enough lattices, this density corresponds 
to the normalized gradient of the lines at fixed 
l3 at the lowest zeros. Hence from the measure- 
ments on a G'* lattice we estimate the transition 
to occur between /3c ~ 0.185 and 0.19. 

On larger lattices, we observe a change in the 
distribution of the zeros. 

For l3 = 0.195 on the 8'* lattice and for /3 < 
0.200 on the lO'* lattice, the lowest zeros are imag- 
inary (at l3 = 0.195 on the lO'* lattice some of the 
lowest zeros fiuctuate very slightly away from the 
imaginary axis but this could be due to inade- 
quate ensemble sampling). At /3 = 0.210 on the 
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Figure 1. The lowest zeros on the G'* lattice plotted 
m sequence. The crosses correspond to [3 = 0.160, 
crossed circles to [3 = 0.185, pluses to [3 = 0.190, 
asterisks to [3 = 0.195, open circles to [3 = 0.200, 
filled diamonds to [3 = 0.205, double circles to 
[3 = 0.210, open diamonds to [3 = 0.215 and the 
crosses on the right hand side to [3 = 0.220, 

lO'* lattice, only the lowest of the small zeros is 
imaginary and at /3 = 0.220 all the zeros have 
gained a small real part. Note that only the zeros 
close to the origin are determined with relatively 
small error. 

Because we did not observe this change in the 
distribution within the above range of [3 on the 
G'* lattice, we performed a simulation at /3 = 0.25 
on a G'* lattice. After a couple of hundred mea- 
surements the lowest zero obtained a small real 
part, but it became complex again after ~ 1000 
measurements. 

Figure 2 shows the lowest zeros on the lO'* 
lattice plotted in sequence for [3 < 0.21. At 
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Figure 2. The lowest zeros on the lO'* lattice plot- 
ted m sequence for [3 < 0.21, The crossed circles 
correspond to [3 = 0.185, the asterisks to [3 = 
0.195 and the open circles to [3 = 0.200. 

j3 = 0.185 the transition clearly looks first or- 
der, hi [3 = 0.195 it appears to be close to first 
order. At that [3 we investigated the scaling be- 
haviours of the lowest zero and of the separations 
between the first and second, and the second and 
third zero. The separations scale with respect to 
the 8^* and lO'* lattices as V" with a just greater 
than —1, again consistent with a behaviour which 
is slightly weaker than first order. The lowest 
zero on the other hand scales with a < — 1. At 
j3 = 0.2, the data from the lO'* lattice indicate 
that p{Q) = 0. 

However, in view of the large finite size effects, 
these conclusions may change when larger vol- 
umes are simulated. 

It is clearly important to study this change in 
the distribution of the zeros on larger lattices. 



The migration of the zeros from the imaginary 
axis could be a signal for a change in the critical 
behaviour of the system. On the other hand it 
could be a result of finite size effects, i.e. our lat- 
tice is too small for the onset of finite size scaling. 
If not, it could be that via finite size scaling the 
real parts of the lowest zeros scale to the imagi- 
nary axis if the critical mass is zero. 

4. Conclusions 

The fact that there is a tricritical point [4] ap- 
pears to be refiected in the rapid increase in mag- 
nitude of the lowest zeros in that region and in 
the change in their distribution. We have evi- 
dence that finite size effects are large in the weak 
coupling sector. 

We intend to perform simulations on larger lat- 
tices to investigate these points in greater detail. 
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